Typeset with jpsjS.cIs <ver.l.O> 



Full Paper 



Cooperative Effect of Coulomb Interaction and Electron-Phonon Coupling on the 
Heavy Fermion State in the Two-Orbital Periodic Anderson Model 



Keisuke Mitsumoto^ * and Yoshiaki Ono^'^ 

^ Graduate School of Science and Technology, Niigata University, Niigata 950-2181, Japan 
^ Department of Physics, Niigata University, Niigata 950-2181, Japan 
^ Center for Trans disciplinary Research, Niigata University, Niigata 950-2181, Japan 



o 

(N 



We investigate the two-orbital periodic Anderson model, where the local orbital fluctuations 
of /-electrons couple with a two-fold degenerate Jahn- Teller phonon, by using the dynamical 
mean-field theory. It is found that the heavy fermion state caused by the Coulomb interaction 
between /-electrons U is largely enhanced due to the electron-phonon coupling g, in contrast to 
the case with the single-orbital periodic Anderson model where the effects of U and g compete 
to each other. In the heavy fermion state for large U and g, both the orbital and lattice 
fluctuations are enhanced, while the charge (valence) and spin fluctuations are suppressed. In 
the strong coupling regime, a sharp soft phonon mode with a large spectral weight is observed 
for small U, while a broad soft phonon mode with a small spectral weight is observed for large 
U. The cooperative effect of U and g for half- filling with two /-electrons per atom n/ = 2 is 
more pronounced than that for quarter- filling with Uf — 1. 

KEYWORDS: heavy fermion, periodic Anderson model, orbital degeneracy, Coulomb interaction, electron- 
phonon coupling, Jahn-Teller phonon, dynamical mean-field theory 
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1. Introduction 

Heavy fermion systems have been intensively investi- 
gated for the last three decades. The effective mass of 
quasiparticles for the systems, estimated from the T- 
linear coefficient of specific heat, the Pauli paramagnetic 
susceptibility and the coefficient of resistivity, is much 
larger than that for simple metals. The strong correla- 
tion effects between /-electrons are considered to play 
crucial roles, where the large entropy due to the local 
magnetic moment of the /-electrons remains down to a 
low temperature To and results in a large effective mass in 
proportion to 1/Tq at low temperature below Tq. When 
the degeneracy is resolved by external magnetic fields, 
the effective mass is expected to be suppressed. Such a 
magnetic suppression of the effective mass is widely ob- 
served in heavy fermion compounds, justifying the heavy 
fermion state of magnetic origin. 

Recently, another class of heavy fermion systems has 
been observed in filled skutterudite compounds RT4X12 
with R=rare earth, T=Fe, Ru, Os, and X=pnictogen, 
where the R ion is surrounded by an oversized cage 
made of X and T atoms. SmOs4Sbi2 shows a large spe- 
cific heat coefficient 7 = 820mJ/K^mol which is al- 
most independent of applied magnetic fields, suggest- 
ing a heavy fermion state of nonmagnetic origin such as 
charge, valence and phonon degrees of freedom. Sm ion 
in SmOs4Sbi2 shows mixed valence state of Sm^"*" and 
Sm^''" as observed in X-ray spectroscopy experiments.^' 
In the Sb-NQR experiment, a peak in I/T2 has been at- 
tributed to some electrical fluctuations, probably local 
charge fluctuations.^-* The specific heat and the ultra- 
sonic measurements have revealed that the crystalline- 
electric-field (CEF) ground state of Sm'^^ ion is Tg? quar- 
tet in T/i site symmetry.^' However, another specific 
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heat measurement claimed the Ts doublet CEF ground 
state and there is still controversy about the CEF state. 

PrOs4Sbi2 with Fi singlet ground state also shows 
a large specific heat coefficient 7 = 750mJ/K^mol to- 
gether with a large jump in the specific heat AC/Tc ^ 
500mJ/K^mol at the superconducting transition temper- 
ature Tc — 1.85K.^^ In the ultrasonic measurements, re- 
markable frequency dependence of the elastic constant 
(ultrasonic dispersion) around 30K has been observed in 
PrOs4Sbi2 and has been attributed to large amplitude 
local vibrations (rattling) of the Pr ion in the cage.^^ 
In addition, PrOs4Sbi2 shows an anomalous softening of 
the elastic constant below lOK down to Tc.^^ The soft- 
ening is well accounted for by quadrupolar susceptibility 

(2) 

due to /-electrons in the Fi — F^ ^ CEF state except for 
an extra softening below 3K,^^ where the coupling be- 
tween the quadrupolar fluctuations and the rattling may 
play important roles for the extra softening and also for 
the heavy fermion behavior down to Tc. Similar proper- 
ties of the ultrasonic dispersion and/or the anomalous 
softening have been observed also in other filled skut- 
terudite compounds such as ROs4Sbi2 (R=La, Sm, Nd) 
and RFe4Sbi2 (R=La, Pr).5.i"-i3) 

After the pioneering works of Holstein,^**^ the prob- 
lem of coupling between local lattice vibrations (Einstein 
phonons) and local charge fiuctuations of conduction 
electrons has been investigated by many authors. Yu and 
Anderson originally predicted that the strong electron- 
phonon coupling causes an effective double- well potential 
for oscillating ions which are responsible for the anoma- 
lously large resistivity and Debye- Waller factor observed 
in A15 compounds such as VaSi and NbaGe.-^^' As a 
strong coupling fixed point, the two- level Kondo systems 
were investigated to describe a heavy-fermion like be- 
havior observed in such compounds. Recently, the 
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precise studies for the electron-phonon impurity mod- 
gjgi9-22) g^j^j ^YiQ multi-level impurity models^"^' have 
been done by using the numerical renormalization group 
(NRG) method and have confirmed the previous re- 
sults including the double-well potential formation in the 
strong coupling regime. The similar approaches have also 
applied to solve the lattice models such as the Holstcin 
modeP^-* and the Hubbard-Holstein modeP^-* in infinite 
dimensions on the basis of the dynamical mean-field the- 
ory^''^ (DMFT) and have revealed that the system be- 
comes bipolaronic insulator when the electron-phonon 
coupling is larger than a critical value, g > gc, where 
the effective potential becomes double- well type. More 
recently, Hattori and Miyake have discussed the ultra- 
sonic dispersion in the Holstein model by using the self- 
consistent ladder approximation. 

To discuss the effect of local lattice vibrations on 
the rare-earth and actinide compounds, the periodic 
Anderson-Holstein model, whore local phonons couple 
with valence fluctuations of /-electrons which interact 
with each other via the Coulomb interaction and hy- 
bridize with conduction electrons, was extensively stud- 
ied in early theories for mixed-valence systems such as 
Sm chalcogenides.^^) As the strong correlation effect due 
to both the Coulomb interaction and the electron-phonon 
coupling is crucial for describing the heavy-fermion state 
in these systems, we need reliable and nonperturbative 
approaches such as the DMFT. Recently, we have stud- 
ied the periodic Anderson-Holstein model by using the 
DMFT.30-32) What we found are as follows: (1) In the 
strong electron-phonon coupling regime, g ^ gc, the sys- 
tem shows an anomalous heavy fermion behavior which 
is accompanied by a large lattice fluctuation and an ex- 
treme phonon softening. (2) A simple harmonic potential 
for ions for g ^ gc changes into an effective double-well 
potential for g ^ gc- (3) The pairing interaction between 
the conduction electrons has a maximum at g « gc- (4) 
The heavy fermion state due to the electron-phonon cou- 
pling is realized in the wide range of the f-electron num- 
ber n/, while that due to the Coulomb interaction is real- 
ized in the narrow range near the half- filling n/ ~ 1. (5) 
The effect of the electron-phonon coupling on the heavy 
fermion state and that of the Coulomb interaction are 
compete with each other. 

In the rare-earth and actinide compounds, orbital fluc- 
tuations of /-electrons are considered to play crucial 
roles in determining the heavy-fermion behavior. '^^^ Re- 
cently, Hotta studied the multiorbital impurity Ander- 
son model in which the orbital fluctuations of /-electrons 
couple with Jahn- Teller (JT) phonons by using the NRG 
method and observed the remarkable quasi-Kondo be- 
havior due to the dynamical JT effect.''"') As for the 
lattice model, we have studied the two-orbital periodic 
Anderson model coupled with JT phonons by using the 
DMFT and have found that the heavy fermion state is 
realized due to the cooperative effect of the Coulomb in- 
teraction and the electron-phonon coupling. This is a 
striking contrast to the case with the periodic Anderson- 
Holstein model where the both effects are compete with 
each other as mentioned above. ^°~^^) However, the effects 
of the degeneracy of JT phonon mode and the /-electron 



filling have not been discussed there. 

In the present paper, we investigate the two-orbital 
periodic Anderson model, where JT phonons with a 
two-fold degenerate mode couple with orbital fluctua- 
tions of /-electrons which interact with each other via 
the Coulomb interaction and hybridize with conduction 
electrons, by using the DMFT to elucidate the coopera- 
tive effect of the Coulomb interaction and the electron- 
phonon coupling on the heavy-fermion state in the cases 
with half- filling with two /-electrons per atom n/ = 2 
and quarter-filling with nf = 1. 

2. Formulation 

Our model Hamiltonian is given by, 

H = ^eucli^cua + ef^fli^fiia 



with 



ila 

+ ^ X! "-/^^t^/iU + ^' nfiiahfi2a' 
il iaa' 

+ ^ filafi2a' fi'^-cr' fi2a 
iaa' 

+ 91 ^{bli + bii)Ti^ + 52 ^{ba + bi2)fi- 

i i 

iv 



(/ilcr/ilo' //2cr/j2cr): 



(1) 



(2) 



(3) 



where cji^ ifua) ^® ^ creation operator for a conduction 
(c-)electron (/-electron) with orbital l{— 1,2) and spin 

(T(=t,i) at site i, and fifua = fl^fua- is a creation 
operator for a Jahn- Teller phonon with mode v{= 1,2) at 
site i, where the normal coordinates for the Jahn- Teller 
phonons are given by Q^^ = I / ^/2uJ^{hil, + b\^). ek, £/, 
and V are the dispersion of c-electron, the atomic level 
of /-electron, and the c-/ hybridization, respectively. In 
the model eq.(l), the Jahn- Teller phonons with the fre- 
quency ujQv couple with the local orbital fluctuations of 
/-electrons, Tix and Tiz, via the electron-phonon coupling 
g^-^^^ The model eq.(l) also includes the Coulomb inter- 
action between /-electrons: the intra and inter orbital di- 
rect Coulomb U and U' and the exchange coupling J. For 
simplicity, we assume wqi = Wo2(= wq), gi = g2{= g), 
U = U' and J = in the present paper. 

To solve the model eq.(l), we use the DMFT^'^) in 
which the model is mapped onto an effective single im- 
purity two-orbital Anderson model coupled with Jahn- 
Teller phonons. The local Green's function G/io-(iw„) 
and the local self-energy T.ia{iujn) for the /-electron sat- 
isfy the following self-consistency conditions: 
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where p{€) is the density of states (DOS) for the c- 
electron, p(e) = 5Zk'^(^ ~ ^k)- In the above equation, 
Gfi{iLOn) is the bare Green's function for the effective 
impurity Anderson model with U = g = in a.n effec- 
tive medium which will be determined self-consistently. 
The effective impurity Anderson model with finite U 
and/or g is solved by using the ED method for a finite- 
size cluster to obtain G//cr(iw„) and E;cr(i^n) at T = 
Q_26, 27, 30-32, 35-37) ^j^g present study, we use 4 site clus- 
ter and the cutoff of phonon number is set to be 9 for 
each Jahn- Teller mode. We assume a semielliptic DOS 
with the bandwidth W = 1, p{e) — — /tt, and we 
set V = 0.1 and wo = 0.01. We concentrate our attention 
on the nonmagnetic state with ^/^(ia'n) = S(ia;„) at half- 
filling with Uf = 2 and quarter- filling with ri/ = 1, where 
nf = '}2i„{nfiia-) is the average number of /-electrons per 
atom. 

3. Results 

3.1 Renormalization factor 

In Fig.l, we plot the renormalization factor Z = 
— as a function of g for several values 

of U at half- and quarter-filling. When g = 0, the ef- 
fective mass of the quasiparticle m* /m = increases 
with increasing U and then the heavy fermion state with 
m* /m ^ 1 is realized in the strong correlation regime. 
When the electron-phonon coupling increases, the effec- 
tive mass increases for all values of U resulting in the 
heavy fermion state due to the cooperative effect of the 
Coulomb interaction and the electron-phonon coupling. 
This is a striking contrast to the case with the single- 
orbital periodic Anderson- Holstcin model where the both 
effects are compete with each other as mentioned be- 
fore.3°-32) The effective mass m* / m for half-filling is con- 
siderably larger than that for quarter-filling as shown in 
Figs.l (a) and (b). This is due to the effect of the local or- 
bital fluctuation which is largely enhanced for half-filling 
as compared to quarter-filling as mentioned later. 

3.2 Local orbital fluctuation 

Fig.2 shows the local orbital fiuctuation written by 

(T|) = (T^) = ((n/a-%i2)') (4) 

with fifii = nfii^ as a function of g for several values 
of U, where (r^) = {tx) = in the absence of orbital 
order and (t|) = (rj) due to the symmetry of eu repre- 
sentation. When .9 = 0, (r^) increases with increasing U 
as previously obtained in the two-orbital periodic Ander- 
son model. When g increases for U ^ 0, (t|) steeply 
increases at half-filling, while it gradually increases at 
quarter-filling except for U = as shown in Figs. 2 (a) 
and (b). Remarkably, the increase in the local orbital 
fluctuation in the strong coupling regime is largely en- 
hanced due to the Coulomb interaction for half-fllling in 
contrast to the case for quarter-filling where it is sup- 
pressed. The difference between half-fllling and quarter- 
fllling in the local orbital fluctuation is due to that in the 
double-occupancy probabilities as shown below. 




g 

Fig. 1. The renormalization fax;tor Z as a function of electron- 
phonon coupling constant g for several values of the Coulomb 
interaction U at half-filling (a) and quarter-filling (b). 




g 

Fig. 2. The local orbital fluctuation (r|) as a function of g for 
several values of U at half- filling (a) and quarter- filling (b). 

3.3 Double occupancy probabilities 

The local orbital fluctuation (t|) given in eq.(4) can 
be rewritten by 

i^z) = {"nfr) + 2 ^{nfil^hfUi) - 2 ^{nf^lanf^2c^') (5) 

I <7<7' 

with fifi = "^infii. To discuss the local orbital fiuc- 
tuation in detail, we also plot the double-occupancy 
probabilities as functions of g for several values of U 
at half- and quarter-filling in Fig. 3, where {fifiianfas) , 
{''T'fiicrfifiia-) {nfiiaf^fQ^) are on-site correlation func- 
tions for intra-orbital spin-singlet, inter-orbital spin- 
singlet and inter-orbital spin-triplet, respectively. When 
(? = 0, the double- occupancy probabilities coincide to 
each other due to the spin-orbital symmetry with U = 
U', and then, they are suppressed due to the electron cor- 
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Fig. 3. (Color online) The double-occupancy probabilities, 
(nfiicrhfiii^) (solid lines), (hfn^nfii^) (dashed lines) and 
{''^ filer'"' fila-) (dotted lines), as functions of g for (7 = and 0.6 
at half-filling (a) and quarter-filling (b), respectively. 

relation effect, especially in the case with quarter-filling 
as shown in Fig. 3 (b). On the other hand, when g ^ 0, 
the degeneracy is resolved, and then, the intra-orbital 
spin-singlet correlation {fifuanfas) is enhanced, while 
the inter-orbital spin-triplet correlation {hfUfjhjjj^) and 
the inter-orbital spin-singlet correlation {nfuanfifa) ^-re 
suppressed due to the effect of g. As the results, (r^) 
given in eq.(5) is largely enhanced for half-filling with 
large values of the double-occupancy probabilities, while 
it is slightly enhanced for quarter-filling with small values 
of those as shown in Fig. 2. 

3.4 Local charge fluctuation 

The local charge fluctuation of /-electrons, i. e., the 
valence fluctuation is written by 

I 

+ 2^{hfii„hfi2a'), (6) 

and is plotted as a function of g for several values of 
U at half- and quarter-filling in Fig. 4. When g ^ 0, 
the valence fiuctuation decreases with increasing U to- 
gether with decreasing double-occupancy probabilities 
for both half- and quarter-filling due to the electron cor- 
relation effect. When g ^ 0, the valence fluctuation is 
suppressed or almost constant due to the effect of g in 
contrast to the case with the single-orbital periodic An- 
derson model where the valence fiuctuation is coupled 
with local phonons and is enhanced due to the effect of 
^ 30-32) note that the sign of the final term (the inter- 
orbital correlation) in eq.(6) is opposite to that in eq.(5), 
and then, the changes in the intra- and inter- orbital cor- 
relations cancel with each other (see Fig. 3) resulting in 
the decrease in the valence fluctuation with g in contrast 
to the increase in the local orbital fluctuation. 




(b) quarter-filling 




g 



Fig. 4. The valence fluctuation ((n^j — (ri/i))^) as a function of 
g for several values of U at half-filling (a) and quarter-filling (b) . 




°0 0.01 0.02 

g 

Fig. 5. The local moment (S'^) as a function of electron-phonon 
coupling constant g for several values of U at half-filling (a) and 
quarter- filling (b). 

3.5 Local moment 

The local moment, (S^) = (Sf) = {Si + Sfy + Sf,), 
written by 

is plotted as a function of g for several values of U at half- 
and quarter-filling in Fig. 5. When g — 0, (S^) is enhanced 
due to the effect of U. When g increases, (S^) decreases 
with decreasing (increasing) spin-triplet (singlet) corre- 
lation as shown in Fig. 3. The decrease in (S^) with g for 
half-filling is larger than that for quarter-filling, as the 
double-occupancy probabilities for half-filling are larger 
than those for quarter- filling (see Fig. 3). 
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Fig. 6. The normalized local lattice fluctuation (Q^) / {Q'^)o, as 
a function of g for several values of U at half-filling (a) and 
quarter-filling (b). 



3.6 Local lattice fluctuation 

The local lattice fluctuation is defined by {Ql) = 
{{Q^. - {Q^.)?) and {QD = {Ql){= (Q^)) due to the 
symmetry of Eg modes. In Fig. 6, we plot the normal- 
ized local lattice fluctuation, {Q^)/{Q^)o, where {Q^)o = 
l/2ijjQ is the value for the zero-point oscillation with 
£1 = 0. When g increases, (Q^) monotonically increases 
for all values of U. Remarkably, the increase in (Q^) 
observed in the strong coupling regime is largely en- 
hanced due to the Coulomb interaction especially for 
half-filling, where the orbital fluctuation coupled with 
the local phonons is largely enhanced. This is a strik- 
ing contrast to the case with the single-orbital periodic 
Anderson model where (Q^) is suppressed due to the 
Coulomb interaction as well as the valence fluctuation as 
mentioned before. 

3.7 F- electron spectral function 

Fig. 7 shows the /-electron spectral function as a func- 
tion of uj at half-filling, where we also plot the analytic 
DOS for [/ = g = in Fig.7 (a) and the DOS eval- 
uated by a Lorentzian broadening with width of 0.1|a;| 
for finite U and/or g in Figs. 7 (b)-(d). For U = g = 0, 
the Fermi level is in the c-f hybridization gap A^f as 
shown in Fig.7 (a). For U — 0.6 and g = shown in 
Fig.7 (c), the c-f hybridization gap is renormalized as 
^ ZAcf with the renormalization factor Z ^ 0.099 due 
to the strong electron correlation effect, i. e., the system 
becomes Kondo insulator. In this case, we also observe 
the lower and upper Hubbard bands around uj — ±U /2. 

For t/ = and g = 0.018 shown in Fig.7 (b), the 
c-f hybridization gap is renormalized as ~ ZAcf with 
Z ~ 0.048 due to the effect of strong electron-phonon 
coupling. We note that, in the case with large phonon 
frequency, ujq ^ W, the system can be described by 
the effective negative-J/cg Anderson model with \Ucs\ — 
2g'^/u!o, where the lower and upper bands are formed 
around ±|[/off|/2. In the present case with small wo ^ 
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Fig. 7. (Color online) The f-elcctron spectral function (vertical 
bars) and the DOS (dotted lines) as functions of ui for several 
values of U and g at half-filling, where the Fermi level is set to 
zero. 



W, however, the clear lower and upper Hubbard bands 
are not observed, instead incoherent sub-bands are ob- 
served together with the quasiparticle bands. For U = 0.6 
and g = 0.004 shown in Fig.7 (d), the c-f hybridiza- 
tion gap becomes narrower (~ ZAcf with Z ~ 0.065) 
than that for U = 0.6 and g = (see Fig.7 (c)) due 
to the cooperative effect of the Coulomb interaction and 
the electron-phonon coupling. In this case, the high en- 
ergy structure of the lower and upper Hubbard bands 
observed for U = 0.6 and g = is almost unchanged by 
the effect of g. 

In Fig. 8, we plot the /-electron spectral function as 
a function of uj at quarter half-filling for U ^ g — 
(Z = 1), C/ = and g = 0.021 [Z - 0.20), U = 0.6 
and g = (Z - 0.11), and U = 0.6 and g = 0.008 
{Z ^ 0.066). The specific features of the /-electron spec- 
tral function and the DOS for quarter-filling shown in 
Fig. 8 are almost the same as those for half-filling shown 
in Fig.7 except that the DOS is asymmetric in contrast 
to the case with half-filling where the DOS is symmet- 
ric. Then, the system becomes heavy fermion metal with 
large effective mass m* /m — Z^^, instead of the Kondo 
insulator with narrow c-f hybridization gap ^ ZAcf re- 
alized for half-filling. 

3.8 Phonon spectral function 

Finally, we calculate the local phonon Green's func- 
tion defined by D^{t) = — {Trbi„{T)bf^{0)) , whose spec- 
tral function is plotted as a function of lu for half-filling in 
Fig. 9 and for quarter-filling in Fig. 10, respectively. When 
g increases, the lowest excited energy shifts to low en- 
ergy as previously observed for the single-orbital peri- 
odic Anderson-Holstein model. ■^""^^^ For a small value of 
U = 0.2, the lowest excited state has a dominant contri- 
bution to the spectral weight. On the other hand, for 
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Fig. 9. The phonon spectral function as a function of ui at half- 
filling for several values of g with U = 0.2 (a) and U = 0.6 (b). 



Fig. 8. (Color online) The f-electron spectral function (vertical 
bars) and the DOS (dotted lines) as functions of ui for several 
values of U and g at quarter-filling, where the Fermi level is set 
to zero. 



a large value oi U = 0.6, the spectral weight of the 
lowest excited state is relatively small and a consider- 
able amount of the spectral weight remains around the 
original phonon level. This result is consistent with that 
obtained for the Hubbard-Holstein model in the strong 
coupling regime, where the sharp soft phonon mode with 
a large spectral weight is observed for small- 1/, while the 
broad soft phonon mode with a small spectral weight is 
observed for large-J/.^®-' 

To see the energy shift in more detail, we show the 
lowest excited energy ujq in the phonon spectral function 
as a function of g in Fig.ll. Remarkably, the softening, 
i. e., the decrease in wq, observed in the strong coupling 
regime is largely enhanced due to the Coulomb inter- 
action especially for half-filling in contrast to the case 
with the single-orbital periodic Anderson- Holstein model 
where the softening is suppressed due to the Coulomb in- 
teraction. ^""^^^ 

4. Summary and Discussions 

In Summary, we have investigated the two-orbital pe- 
riodic Anderson model coupled with the two-fold degen- 
erate Jahn- Teller phonon by using the dynamical mean- 
field theory, and have found that the heavy fermion state 
of nonmagnetic origin is realized due to the cooperative 
effect of the Coulomb interaction U and the electron- 
phonon coupling g. The specific features of the heavy 
fermion state for large U and g are as follows: (1) The 
local orbital and lattice fluctuations are enhanced, while 
the local charge (valence) and spin fluctuations are sup- 
pressed. (2) The sharp soft phonon mode with a large 
spectral weight is observed for small U, while the broad 
soft phonon mode with a small spectral weight is ob- 
served for large U. (3) The cooperative effect for half- 
filling with nj = 2 is more pronounced than that for 
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Fig. 10. The phonon spectral function as a function of at 
quarter-filling for several values of g with U = 0.2 (a) and 
U = 0.6 (b). 



quarter- filling with n/ = 1. 

In the absence of g, when U increases, the effective 
mass increases with increasing spin and orbital fiuctua- 
tions together with decreasing charge fluctuation. Then, 
the heavy fermion state of magnetic origin with large 
(small) spin and orbital (charge) fluctuations is realized 
for large U due to the strong correlation effect. When we 
increase g for a flxed value of U, the effective mass further 
increases with increasing orbital and lattice fluctuations 
together with decreasing spin and charge fluctuations. 
Then, the heavy fermion state of nonmagnetic origin with 
large (small) orbital and lattice (spin and charge) fluctu- 
ations is realized for large U and g due to the cooperative 
effect of the strong correlation and the strong coupling. 
In this heavy fermion state, the phonon spectral function 
shows the broad soft phonon mode while the / electron 
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Fig. 11. The lowest excited energy uio in the phonon spectral 
function as a function of g for several values of U at half-filling 
(a) and quarter- filling (b). 



spectral function is almost unchanged due to the effect 
of g except further narrowing of the low energy quasi- 
particle bands. 

In the two-orbital periodic Anderson model with the 
coupling g between the local orbital fluctuation and the 
Jahn- Teller phonon, the orbital fluctuation is enhanced 
due to the both effects of U and g, and then the heavy- 
fermion state is realized due to the cooperative effect. 
This is a striking contrast to the case with the single- 
orbital periodic Anderson model with the coupling g be- 
tween the local charge fluctuation and the local phonon, 
where the spin (charge) fluctuation is enhanced (sup- 
pressed) due to the effect of U, while the charge (spin) 
fluctuation is enhanced (suppressed) due; to the effect 
of g, and then the effects of U and g on the heavy- 
fermion state compete with each other. ^""^^^ As the ab- 
solute value of the local orbital fluctuation for = 2 
is larger than that for nf = 1, the cooperative effect of 
U and 5 for n/ = 2 is more pronounced than that for 
Uf = 1. 

Although the present model is the simplest version, the 
multi orbital periodic Anderson model coupled with the 
Jahn- Teller phonons has the potential to account for the 
experimental observations. The heavy fermion state of 
nonmagnetic origin due to the cooperative effect of U and 
g seems to be consistent with the magnetically robust 
heavy fermion state observed in the filled skutterudite 
SmOs4Sbi2. The enhancement of the orbital fluctuation 
and the anomaloiis phonon softening due to the cooper- 
ative effect of U and g seemes to be responsible for the 
extra softening of the low temperature elastic constant 
which is not reproduced only due to the quadrupolar sus- 
ceptibility in the CEF state observed in PrOs4Sbi2. To 
discuss these compounds in detail, however, we need to 
consider more realistic models including the effects such 
as the spin-orbit interaction, the Hund's rule coupling 
and the crystalline-electric-fleld. Within the j-j coupling 



scheme, Hotta has developed the realistic three orbital 
Anderson model with r7 doublet and Fg quartet CEF 
states simulating filled skutterudite compounds. It is 
interesting to extend our approach to such realistic mod- 
els and is now under the way. 
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